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Abstract: We introduce Q-space, the tensor product of an index space with a pri¬ 
mary space, to achieve a more general mathematical description of correlations in terms 
of (/-tuples. Topics discussed include the decomposition of Q-space into a sum-variable 
(location) subspace § plus an orthogonal difference-variable subspace B, and a systema¬ 
tisation of g-tuple size estimation in terms of p-norms. The “GHP sum” prescription for 
g-tuple size emerges naturally as the 2-norm of difference-space vectors. Maximum- and 
minimum-size prescriptions are found to be special cases of a continuum of p-sizes. 


1 Correlations in P-space 

Traditionally, particles emitted by high-energy hadronic or heavy ion collisions have been 
visualised in terms of a collection of points populating what we call primary space P or 
P-space, with each particle i represented by a d-dimensional vector x* = (x t \ , x* 2 - • • • ,Xid)- 
Examples of P-spaces are three-momentum, with x,; = ( Pi x ,Piy,Piz ) and rapidity-azimuth, 
with Xj = (yi,<j>i). Particle correlations can be studied either by binning P with a suit¬ 
able partition (“coarse-graining”), or by analysing distributions of relative distances between 
primary vectors x, directly. In this contribution, we focus exclusively on the latter approach. 

Correlations are a matter of definition; typically they are deviations of joint g-particle distri¬ 
butions from a pre-defined null hypothesis or reference process[ 1] such as a uniform distribu¬ 
tion or a (/-fold convolution of the differential one-particle distribution. Under a “dilute-fluid” 
assumption that correlation strength decreases with particle cluster ((/-tuple) size, low-order 
(/-tuples — particle pairs, triplets, quartets etc. — are usually studied as a first approxima¬ 
tion, with higher-order (/-tuples as perturbations. In general, one selects out of N particles all 
possible combinations 2 of (/-tuples [q = 1,2,3,...) for statistical analysis, using for example 
cumulants[ 2] as differential correlation measures. The characterisation of (/-tuples is therefore 
fundamental to correlation analysis. 

The simplest properties of a given (/-tuple are its location and size. Location can be defined, 
for example, as the centre of mass of the q particles, 

x = - X Xi > W 

or as the location of any one of the particles. 

In its simplest incarnation, the mathematical realisation of size should be a nonnegative 
real number which reduces to zero whenever all q particles occupy the same position in 

1 in: Correlations and Fluctuations in QCD, Proceedings of the 10th International Workshop on Multi¬ 
particle Production, Crete, June 2002, edited by N.G. Antoniou, F.K. Diakonos and C.N. Ktorides, World 
Scientific (2003), pp. 386-403. 

2 In analytical manipulations, it is more convenient to consider all N\/(N — q)\ ordered g-tuples rather than 
the N\/(N — q)\q\ unordered ones. 
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P: in other words, size should be a norm based on relative coordinates. Contrary to naive 
expectation, the best prescription for this is not immediately obvious. While a 2-tuple’s size 
is clearly described in terms of the distance |x, — x y -|, (/-tuples of higher order permit a range of 
choices. In Ref.[3] this problem was addressed in terms of different “topologies”, summarised 
pictorially in Fig. 1 for some representative 4-tuples. 



Figure 1: An event, shown as a collection of N points in the P-space (y,4>). Also shown are 
examples of three topologies used to quantify 4-tuple size. 

For every topology, interpair distances can be combined in several ways to form different 
size estimators. The q(q — l)/2 interpair distances making up the GHP topology can, for 
example, be summed to yield the GHP sum size; or size can be based on the largest of these 
distances, resulting in the GHP max size estimate. Particle number within d-dimensional 
spheres centered on individual particles defines the Star max size, while the Snake Integral 
seeks to quantify size in terms of a linear succession of distances between ordered points. 

Historically, the GHP max and Star max prescriptions were heuristic inventions within the 
correlation integral literature: the q= 2 simplest case[4] was extended in Refs.[5] and[6] to the 
GHP max prescription, 

C (Gm> \i) = {No. q-tuples (R, • • •, i q ) with all |x in - x im | < 1} , (2) 

and in Refs. [7, 8, 9] to the Star max, 

i =i y =1 

(with 0[x) the Heaviside function), which can, on multiplying out the (q— 1) inner sums and 
using 9(£— |xj—Xj|) = 9(1— maxj(|xj—Xj|)), be made to exhibit the Star max prescription 
on the (/-tuple (i,ji, ■ ■ ■ ■ j q - \) explicitly, 

1 N N 

cf tarM) = Y e ^~ m i x * - x a i}) • ( 4 ) 

*=i ji,-j 9 -i=i 

Factorial extensions were published in Ref. [3]. Note that, in summing indiscriminately over 
all g-tuples, the above correlation integrals implicitly assume that correlation structure is 
independent of location. 

2 Q-space by example 

The concept of Q-spaces is not new: In high-energy hadronic collisions, two-particle corre¬ 
lations have long been visualised in two-particle spaces[10]. Figure 2 illustrates by means 





of a simple example how a Q-space is constructed from P-space. A typical event in a one¬ 
dimensional {d= 1 ) primary space is represented by the dots on the lines below the xi axis 
and to the left of the X 2 axis. Particle pairs are then represented by all possible dots in the 
< 5=2 space as shown by representative dashed lines . 3 

Each pair in this example is represented by a vector X^ as shown. The location corresponds 
to the component vector X 5 along the diagonal, and the 2 -tuple size to the magnitude of X^ 1 , 
since |X s | = \x\+X 2 \/\/2 = v/SxandlX 15 ! = \x\ — ,T 2 |/\/ 2 . This can be derived algebraically 
by using the Q-space representation X'? = aqui + X 2 U 2 in terms of the index unit vectors 
Ui 2 , which are rotated to the basis vectors hg 2 shown in Fig. 2 to find X ,s = hi (hi • X®) 
and X D = h 2 (h 2 -X Q ). 



Figure 2: Example of Q-space for q=2, d=l. The primary space event, represented on the 
margins, maps onto particle pairs in the Q-space square. Vector X^, representing a particular 
pair, decomposes into a location vector X ,s and a size vector X D . Unit vectors U 1.2 are rotated 
to hi and h 2 which span sum space and difference space respectively. 

Correspondingly, 3-tuples can be visualised for one-dimensional P-spaces as shown in Fig. 3. 
Here, one particular ordered triplet is shown in all its exchange-permutation incarnations, 
all of which represent the same unordered triplet; the associated symmetry is clearly visible 
when viewing the Q-space down the main diagonal as in Fig. 3(b). The set of unit vectors 
(ui, U2, U3) is transformed to 


hi 

= (ui + u 2 + usVVs 

(5) 

h 2 

= (ui — u 2 )/ V2 

( 6 ) 

h 3 

= (ui + u 2 - 2 u 3 )/\/6 

(7) 


where hi points along the main diagonal (shown as the dashed line), while I 12 and I 13 span 
the plane, indicated by the dotted lines in Fig. 3, which is normal to the main diagonal. This 
normal plane exemplifies difference space D 3 , within which only relative coordinates appear, 
while the main diagonal represents the sum space §3 measuring, once again, g-tuple location. 


3 Making up a “pair” from a particle with itself would result in a dot lying on the diagonal line. Such 
usage corresponds to the transition from factorial to ordinary statistics. We ignore associated issues in this 
contribution. 




Figure 3: A 3-tuple (with d= 1) showing (a) a 3-tuple point plus its five permutation- 
symmetric counterparts and the main diagonal (dashed line), and (b) the same system when 
viewed down the main diagonal. Also shown are the corresponding difference vectors X s . 
The dotted lines are coplanar with the difference space D 3 . 


3 Formalism for Q-space 


The general formalism for vectors in Q is now easily understood. The particle vectors x;,z = 
1,... q of a (/-tuple in d-dimensional primary space P are combined into the Q-space-vector 

g 

X« = ^x i u i . (8) 

i =1 


Q is spanned by unit vectors which are the product of unit vectors u= 1,... ,q living in 
index-space I q and the basis vectors of P (implicit in x;). The d-dimensional sum space S q is 
spanned by 





(9) 


and the basis vectors of P. The corresponding index-space projection of X^ onto sum space, 
X s = hi (hi • X^), is given by 


X" = ^5> = ^x, (10) 

v 9 i =1 

where x is the g-tuple CMS as before. Furthermore, the algebraic complement of 8> q is spanned 
by the set of (q — 1) orthonormal vectors 4 in Ig 


1^2 

= (ui — U 2 )/ V 2 , 

(ID 

h3 

= (ui + u 2 - 2u 3 )/\/e 

(12) 

hg 


(13) 


which, together with the P basis, are used to define a basis for difference space D g . The index 
space projection of X^ onto O q is given by 

_ i =2 

4 Any basis set connected to this one by an orthogonal transformation will also do. 


A 


(14) 





Since every £ Q has a unique decomposition into X 5 and X D . Q is the direct sum of § g 
and B g . The relationship between the different spaces can hence be summarised as 

I g ®P = Q = S g ©B g , (15) 

corresponding to the dimensionality relation dim(I g ) x dirn(P) = dim(Q) = dim(S g )+dim(B g ) 
or, explicitly, qxd = qd = d+(q— 1 )d. 

Within this framework, the first and most obvious measure of size is the 2-norm of X D . 
Starting either with Eq. (14) or by subtraction, X^ = ~Xy — X , we find the suggestive form 


<? 

X D = ^ (x, : - x) U,;, 
i —1 


(16) 


which yields 



(17) 


which, apart from the prefactor, is exactly the GHP 2-sum or rms size measure. The GHP 
2 -sum thus appears to be the natural measure of size in Q-space. 


4 Generalised p-sizes 

While the above result represents a strong endorsement of the Q-space approach to mul¬ 
tiparticle correlations, there is more. Generalised p-norms are the extension of Eq. (17) to 
arbitrary p > 1: Given an m-dimensional vector A = (ai,..., a m ), the p-norm of A is defined 
by 

|ai p = [X>ii p ] 1/ 7 (*> 

i =1 

Special cases are the p =2 norm appearing in Eq. (17) and the “max” norm, 

lim ||A|| = max( |cij|). (19) 

p —>oo 1 i 

The condition p > 1, required for ||A|| p to satisfy the Minkowski inequality ||A + B|| p < 
| A| + [B|| p , can be relaxed to p £ R if we do not insist on ||A|| being a norm and allow it to 
be merely a “size measure” or “p-size”. We then also have negative-p sizes and in particular 

lim |A| = min( |a;|), |a.j| / 0 Vi. (20) 

p—> —oo 1 i 

Fig. 4 shows as a simple example for a two-component vector A = (aq, X 2 ) the set of isonorms 
satisfying ||A|| p = 1 = [|xi| p + |a ?2 | p ] 1//p i.e. a set of curves showing which vectors A have 
the same “p-size”. Only positive x\ and aq are shown as ||A|| is reflection-symmetric about 
x\=0 and X2=0. The usual circle for constant 2-norm is complemented by the straight 
diagonal line (or, in the full plane, diamond shape) of p=l and various shapes in between. 
Of particular interest are the max and min size measures shown as the solid line and dashed 
line respectively. 

Based on the above, we can define a “GHP sum p-size” as follows: 

Ixi-xylf 7 ', P# 0 ( 21 ) 

i<j =1 




Figure 4: Example of p-sizes: A vector A = (x\,X 2 ) has many different norms or, more 
generally, p-sizes. Conversely, different sets of vectors A have the same p-sizes, as shown in 
this picture. The usual quarter-circle Pythagorean norm (p= 2) is complemented by various 
convex and concave “isonorm” curves. Of particular interest are the p = +oo (maximum) 
and p = —oo (minimum) sizes. 


which for p > 1 is a norm. Eq. (17) is seen to be the special case ' = yjq ||X' D || 2 , while 

the GHP max definition is the special case p —> oo, 


Sp = max |xj — x,j 
*<l=i 


( 22 ) 


representing the size prescription used in the correlation integral ( 2 ). 


A given Q-space vector A will therefore yield a set of size measures |||A|| p \ p £ M, p ^ 0|. 
which includes the “min”, “max” as well as the usual 2-norms. An infinite set of p-sizes 
of a given (/-tuple is, however, mostly redundant; for practical purposes, a subset such as 
p £ {—oo, 0.5, 2, Too} is probably sufficient. 


The sizes S qtP can be understood in terms of projections in Q-space as follows. Define the 
set of pair plane normal vectors h ij = (u,; — Uj)/y/2; these can be considered to span the 
respective pair s difference space D 2 ’^. Eq. (21) then can be written as 


r « 

c(GHP) r V” 

*<i=i 




i tv 


(23) 


The set of hj/’s are clearly not mutually orthogonal, given that q(q— 1)/2 vectors all live in 
the (q — l)-dimensional difference space indeed, we have h ij ■ h^ = (S^ + 5ji — 5u — 5jk )/2 
and, conversely, Eqs. (11)—( 13) show that the h . q >2 are simple sums of the hjj’s. Fig. 5 shows, 
for the simple case q =3 and ri=l, the normal vectors in the difference space B 3 (equivalent 
to viewing the Q cube of Fig. 3(a) down the main diagonal) as well as the distances from the 
respective pair planes dij = |x* — Xj |/y/2. 


a 





Figure 5: Example (d= 1, < 7 = 3 ) of the nonorthogonal projection ( not decomposition!) of Q- 
space vectors in terms onto “pair plane normals” h,.. ; . The difference space vector X 15 , shown 
as a point in (b), is projected onto the pair plane normal vectors hi 2 , I 123 and I 131 to yield 
the distances dij between ~K D and the respective pair planes |xj — Xj| = 0. The views are 
once again down the main diagonal. 

5 Q-space and other size measures 

The symmetry of Q-space clearly favours the GHP topology over the corresponding Star and 
Snake topologies. It is possible, nevertheless, to accommodate the latter into Q-space. A 
simple ad hoc definition of the Star p-size of the (/-tuple centered on x,; would be 

S( s ">( X .) EE []T X,; - X ,f] 1/P , p + 0. (24) 

3 =1 

This is the p-generalised size measure corresponding to the Star max correlation integrals 
of Eq. (3)—(4). Alternatively, we can, for the Star (/-tuple centered on i= 1, define a vector 
Xgt ar ( x i) = y^xihi an< i from this find the vector 

Xgar(xi) = X Q - Xf tar (xi) = (x 2 - xi)u 2 + • • ■ + (x, - xi)u, ; (25) 

then the size measure (24) can be applied directly in terms of its u j components. Clearly, 
the pair (Xg tar (xi), X^ ar (xi)) is not orthogonal to hi and hence lives neither in § g nor D g ; 
also, this pair will be different for every (/-tuple centre Xj. This is consistent with the fact 
that sizes of Star (/-tuples are, by definition, different for the different centre particles. 

Star and Snake topologies can be accommodated in the Q-space presented here, but, due 
to their obvious asymmetry, do not fit comfortably into the explicitly symmetric Q-space 
framework. Other approaches such as conditional spaces (for the Star) and spaces based 
on strict ordering such as in time series (for the Snake) will probably yield more natural 
interpretations for these cases. 

6 Summary 

It is becoming clear that the concept of Q-spaces is yielding insight into the fundamental 
structure of correlations of point sets. As subspaces of the final-state phase space q = 
N, Q-spaces have a solid theoretical and historical foundation, while offering a structured 
set of numbers characterising (/-tuples, starting with location and size. Generalising simple 
notions of norms to an infinite set of p-sizes, we find that the “max” and 2 -sum sizes are 
just special cases within this wider set. Nonorthogonal “pair planes” projections are found 
to be important, in keeping with the obvious fact that relative distances are also mutually 
constrained and hence dependent. 





Most importantly, the usual 2-norm of the difference space vector is found to be the GHP 
sum size prescription, which therefore should be afforded more attention in higher order 
correlation analyses. 

Extensions such as shape characterisers and higher-order measures of size are easily conjured 
up. 

Acknowledgements 

This work was supported in part by the National Research Foundation of South Africa and 
by the United States Department of Energy. 


References 

[1] H.C. Eggers, in: 30th Int. Symp. on Multiparticle Dynamics, World Scientific (2001) pp. 
291-302; hep-ex/0102005 

[2] A. Stuart and J.K. Ord, Kendall’s Advanced Theory of Statistics , Vol.l, fifth edition, 
Oxford University Press, New York (1987). 

[3] P. Lipa, P. Carruthers, H. C. Eggers and B. Buschbeck, Phys. Lett. B285, 300 (1992); 
H. C. Eggers, P. Lipa, P. Carruthers and B. Buschbeck, Phys. Rev. D48, 2040 (1993). 

[4] P. Grassberger and I. Procaccia, Phys. Rev. Lett. 50, 346 (1983). 

[5] Hentschel and Procaccia, Physica 8D, 435 (1983). 

[6] P. Grassberger, Phys. Lett. A97, 227 (1983). 

[7] G. Paladin and A. Vulpiani, Lett. Nuovo Cimento 41, 82 (1984). 

[8] K. Pawelzik and H. G. Schuster, Phys. Rev. A35,481 (1987). 

[9] H. Atmanspacher, H. Scheingraber and G. Wiedenmann, Phys. Rev. A40, 3954 (1989). 

[10] L. Foa, Phys. Rep. 22, 1 (1975). 



